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Abstract 

We study a model of N complex component symplectic fermions in D spacetime dimensions. It 
has an infra-red stable fixed point in 2 < D < 4 dimensions we refer to as Sp^J . Based primarily 
on the comparison of exponents, we conjecture that the critical exponents for the 3D Wilson-Fisher 

(3) 

fixed point for an O(N) invariant iV-component bosonic field can be computed in the Sp_2N theory. 

(3) 

The 3D Ising model then corresponds to Sp_2- To lowest order, the exponent j3 agrees with known 
results to 1 part in 1000 and is within current error bars. The v exponent agrees to 1% and we 
suggest this is because we only went to 1-loop for this exponent. 
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I. INTRODUCTION 



In 2D there exists a vast variety of well-understood critical points, i.e. fixed points 
of the renormalization group (RG), and for many practical purposes there is a complete 
classification In 3D the known critical points are comparatively rare. The best known 
is the Wilson-Fisher fixed point which can be understood as a linear sigma model for an 
M-component vector field n with the euclidean space action: 

Swf = J d D x Q<V? • d^n + A (n ■ nf^j (1) 

where d^d^ = Y^=i We will refer to the fixed point theory as It describes classical 

temperature phase transitions in a variety of magnetic systems. 



In our recent work[3], motivatedprimarily by the search for a deconfined quantum critical 
point in 2d anti-ferromagnetism 4|, |5| , we considered a symplectic fermion theory with action: 

S x = J d D x (d, X ] d,x + 16vr 2 A \ X ] x\ 2 ) (2) 

where x is an iV-component complex fermionic field, x^X = YliLiXiXi- We refer to this 
model as a symplectic fermion because in terms of real fields, each component can be written 
as x — Vi + ^2, X^ — Vi ~~ iV2 an d the free action is 

S = i J d D x €ij d^id^j (3) 

where e\2 = an d thus has Sp(2) symmetry. The 2N real component theory is invariant 
under T] — > Urj where U is a 2N x 2N matrix satisfying U 1 \nU = where = e ® In- 
Thus the complex A^-component theory has Sp(2N) symmetry. The symplectic fermion is 
known to have important applications in physics, for instance to dense polymers [7|] in 2D, 



and to disordered Dirac fermions in 2 + 1 dimensions [ 

It is well-known that the symplectic fermion is non-unitary. In possible applications to 
2d anti-ferromagnetism, this issue was addressed in j^], so we do not repeat this here since 
the present context is very different. For statistical mechanics, as we will explain below, the 
non-unitarity is of a very different nature. 

As shown in [3J], the (x^x) 2 interactions drive the theory to a new infrared stable fixed 



point, which we will refer to as Sp^ |23|| . In applications to quantum critical spin liquids 



SIB 



in 2d, the x fields, for N = 2, were proposed to describe a deconfined quantum critical point 



of an anti-ferromagnetic spin system where the spin field n is composite in terms of x- 
n = x £?X- In the N=2 component case, a are just the Pauli matrices. As will become clear, 
in the present context there will be no need to generalize the above expression for arbitrary 
N since n will not be a composite field. 

In this article we focus on applications of the critical theories Sp^j to statistical mechanics 
in 3D. The main conjecture is that, for the purpose of calculating the critical exponents, 
the following two models are equivalent 

Off = Sp%, for M > (4) 

We emphasize that we are not claiming a complete equivalence of the above models, but 

(3) 

rather, that some of the M model exponents can be calculated in the symplectic fermion 
theory. This conjecture is based on an argument that Sp 21 y for negative N should correspond 
to M models, however not according to the naive equivalence N = —M/2. (See below 
for explanations.) We then conjecture N = — M based mainly on the agreement of the 
exponents and an additional symmetry argument. Since the Ising model is known to be 

(3) 

in the universality class of 0\ , we compute 3D Ising model exponents, and exponents for 
other models as well, and show that to lowest order they agree exceedingly well with known 
results. 



II. SCALING THEORY AND RG ANALYSIS 

In this section we summarize the results of the RG analysis in [3|]. The beta function, 
computed to 1-loop only, for the N component model in D space-time dimensions is 

^ = (4-D)A + (iV-4)A 2 (5) 

where increasing the length e e corresponds to the flow toward low energies. The above beta 
function has a zero at 

Note that A* changes sign at N = 4. It is not necessarily a problem to have a fixed point 
at negative A since the particles are fermionic: the energy is not unbounded from below 
because of the Fermi sea. Near A* one has that dX/di ~ (D — 4) (A — A*) which implies the 
fixed point is IR stable regardless of the sign of A*, so long as D < 4. 



For some exponents we need a measure of the departure from the critical point; these are 
the parameters that are tuned to the critical point in simulations and experiments: 

S X ^S X + J d D x (mVx) (7) 

Above, m is a mass. The correlation length exponent v is then defined as 

i ~ m- v (8) 

To streamline the discussion, let \X\ denote the scaling dimension of X in energy units, 
including the non- anomalous classical contribution which depends on D. An action S nec- 
essarily has \S\ =0. Using [cPx] = —D, the classical dimensions of couplings and fields are 
determined from [S*] = 0. The exponents are functions of the anomalous dimensions 7 x ,7 m 
of x an d m: 

[ x ] = (£>-2)/2 + 7x , M = l-7m (9) 
The 7 X exponent determines the two point function of the x fields: 

(X f (x)x(0)) ~ — (10) 

Since [£] = —1, this implies v = — = l/[m]. Let us also define the exponent 

/3 = 2[x]/[m]: 

I =1 _ 7m , ^ = ^±^ 

^ 1 - 7m 

The normalization of the above exponents is appropriate for the context of deconfined quan- 
tum criticality, and will be related to conventional normalizations of the exponents of Om 
models in the next section. 

In j^], 7 m was computed to 1-loop and 7 X to two-loops. The computation can be un- 
derstood as essentially the lowest order epsilon expansion around D = 4. We found the 
following x field exponents for any N, D: 

(A-D)(l-N) (4 - D)\l - N) 

7m ~ 2(4 -iV) ' 7x ~ 4(4-JV)2 { ' 

These in turn imply the following: 

2(4 - N) 



\ L — LJ I . \ ■ LJ 1~4 ^ ^ 

9(n - 9V n 2 -4W-H9U ^> 2 (^ - n\ 

P = 



(4 - N)(D(1 - N) + 2(iV + 2)) 
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III. 3D STATISTICAL MECHANICS 



Since we are now mainly interested in the context of statistical mechanics in D spacial 
dimensions, we need to reconsider the definitions of the exponents in order to compare with 

(3) 

what is known about the M models. In analogy with Ising models, consider an action 

S = S crit icai + J d D xt$ t (x) (14) 

where t = (T — T c )/T c , and $ e is an 'energy operator'. The exponent ^ stat is defined by the 
two point function of the spin field a(x): 

and the exponents z/ stat , /3 stat as follows: 

£^t~^ tat , (a)~i Atat (16) 

The spin field is identified as the M-vector n. It is well-known that the above exponents 
satisfy the scaling relation: 

)/2 (17) 

As explained below, the above exponents are related to u, (3 of the previous section by 
frtat = v/2 and /? stat = /3/4. 
'm 



Let Ovp denote the Wilson-Fisher fixed point. First let us point out that there is a 



simple argument that Offi should correspond to Sp_^ . We emphasize that our conjecture 
eq. (j4j) differs by a 2. The naive argument goes as follows. Suppose we try to make 
the functional integral over \ gaussian by introducing an auxiliary field u(x) and the action 
S = Sfree+ / d D x (u X~ A _1 m 2 ). Then the functional integral over , X gives a w-dependent 
functional determinant to the iV-th power. On the other hand, if x were a complex boson, 
one would obtain the same determinant to the (— ]V)-th power. Since the Om model has M 
real bosonic components, then one naively expects the previously mentioned equivalence. 
One reason we have to doubt this naive equivalence is that it is not obvious that one can 
compute correlation functions of the x's from the effective action for u since u is essentially 
the composite field x^X- Below, we will compare the (3 exponent for the Off and Sp^J 
models and show that they do not coincide under the equivalence iV = —M/2. It also seems 
likely that at finite temperature the naive equivalence iV = —M/2 will be spoiled by the 



different boundary conditions (periodic verses anti-periodic) for the x verses n fields which 
are needed to compute the determinant. 

Another, though related, approach to identifying our models at negative N is based on 
Parisi-Sourlas supersymmetry[9|]. However as we now explain this would give the same 
identification as the considerations of the last paragraph. For a recent discussion of this in 
the context of 2D non-linear sigma models see[lo|. In this construction, one considers both 
bosonic and fermionic fields: 

4>= {(p!, - ■ ■ ,4>M+2N,Vir ■ - V2n) (18) 

where (fi a are bosonic, r\{ are fermionic, and both are real. The bilinear <fi ■ (f> — <?Wa + 
j{ e N)ij ViVj nas OSp(M + 2N\2N) symmetry and one can consider a linear sigma model 
with ((f)- (f)) 2 interactions. The Parisi-Sourlas arguments suggest that exponents should only 
depend on M, i.e. are independent of N, at least if this symmetry is preserved. The purely 
bosonic N = model has Om symmetry where <fi ~ n, whereas the purely fermionic model 
M + 2N = has Sp2N symmetry. Chosing M + 2N = 0, one naively has the equivalence 
Om ~ Sp-Mi which again differs by a factor of two from our conjectured equivalence. The 
observations of the last paragraph still apply, and we again point out that our two-loop 
calculations for (3 appear to contradict the supersymmetry argument. 

Since our conjecture Om ~ Sp-2M differs from the above Om ~ S^-m, clearly some new 
and non-obvious arguments are needed in support of it, which we now attempt to provide. 
The main evidence comes from comparison with exponents. As we show below, the lowest 
order exponents agree for Sp-2M agree exceedingly well with the very high order epsilon 
expansion of the Om model. (See the table below.) For instance, for the Ising model, (3 is 
known to be approximately 1.303 ± .006 whereas the 5p_2 model gives (3 = 13/10 = 1.300. 
The Sp-i model on the other hand (N = —1/2) gives (3 = 56/45 = 1.244. 

A symmetry argument supporting our conjecture can be given as follows. First observe 
that the Sp(2N) symmetry of the symplectic fermions has an O(N) subgroup, generated 
by 1 (g) A where A is an N x N antisymmetric matrix. Near D = 2, the fixed points are 
expected to be described by a non-linear sigma model. In terms of the x-fields, the non- 
linear constraint is x^X = !• As before, expressing each component of x i n terms of real 
components T], xi — Vi + "72, X^ = Vi ~ ^V2, etc, then if r] is fermionic satisfying r/ 2 = 0, Vz, 
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the constraint reads 

V1V2 + V2V3 H h V2N-1V2N = 1 (19) 

It is evident that this constraint has the above-mentioned O(N) symmetry. On the other 
hand, if the rj were bosonic, the constraint reads Ylf=i Vi = 1> which has an 0(2N) symmetry. 
Thus, based only these symmetry arguments, the O(N) models could have realizations in the 
symplectic fermion theory Sp-2N- The flip in sign in Sp-2N is then understood as coming 
from the need to flip the statistics of x- 

Another, somewhat different, argument comes from 2D. What we can learn from the 
well- understood 2D case is somewhat limited because the phase structure is quite different 
than D > 2[llJ. Nevertheless, we can give some additional support for our conjecture as 
follows. The O(M) models are only critical for —2 < M < 2 in 2D where they are equivalent 
to loop models. Parameterizing M as M = —2 cos Tig, the Virasoro central charge is known 
to be c = 1 — Q(g — l) 2 /g. The region < g < 1 is referred to a dense loops, and g > 1 as 
dilute. First note that g = 2 in the dilute phase corresponds to M = —2 and has c = —2. 
This value of c is known to correspond to an iV = 1 component symplectic fermion, thus here 
one observes 0_2 ~ Sp2, in accordance with standard Parisi-Sourlas symmetry. However 
there exists another interpretation in the dense phase that coincides with our conjecture. 



The M models in 2D we r e studl ed using supersymmetry in [12, U* " « argued m 
|13l | that the dense loop phase is not generic, i.e. allowing crossings of the loops leads 
to a dangerously irrelevant perturbation which can drive it to a distinct massless phase, 
referred to as a Goldstone phase. Arguments were also given that this Goldstone phase 
is more generic and what one expects in D > 2. The case OSp(l\2), which corresponds 
to 0_i, is a good example since both M,N = 1 are positive, and was studied in detail. 
M = — 1 is realized in the dense loop phase at g = 1/3. The Goldstone phase is expected to 
have massless degrees of freedom described by OSp(l\2)/OSp(0\2) ~ S ' 2 where S ' 2 is the 
supersphere. The non-linear sigma model has the action 

S =ljd 2 x {{d^f + idwdw) (20) 

where is a single component real field and the sigma model constraint is <fi 2 + ^1^2 = 1. 
Using the Grassman nature of the 77's, the constraint can easily be solved as <fi = 1 — \f]\r\2- 
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This leads to the action 



S = j J d 2 x (idriidrj 2 - ^l^cfyid^J (21) 

In D > 2 the interaction is irrelevent and one ends up with a single N = 1 complex 
component symplectic fermion with c = —2. This value of c was verified with the Bethe 
ansatz in 12|. So in this case, 0_i ~ Sp2, which coincides with our conjecture eq. (j4j) 
extended to negative M. 

It should be clear by now that if the explanation of our conjecture is based on the above 
kinds of arguments, then computing exponents for Om models in an Sp_2M model is simply 
a trick and one need not be concerned about the non-unitarity of the symplectic fermion. 

The above arguments are not a complete explanation and the equivalence conjectured in 
the introduction was originally based on the comparison of exponents. Let us identify $ e as 
X*X- Because of the m 2 in eq. (JTj) one finds 

fstat = v/2 (22) 
Supersymmetry implies that \n\ = [x] since they are both in the same multiplet 0, thus 

The exponent r] stat then follows from the scaling equation (1T7T) . which imples 77 sta t = 2j x . 
In 3D the formulas of section II give 

_ 2(4 - N) _ 2N 2 - UN + 33 

V ~ 7-N ' * ~ N^-llN + 2% ^ ' 

Note that eqns. ( I22f23l) in the N — > — oo limit give z/ sta t — > 1 and /? s t a t - > 1/2 which is 
consistent with the M — > oo limit of the Om model and also explains the factor of 4 in eq. 
(El. 



(3) 

The equivalent lowest order epsilon-expansion result in the bosonic M model ispj 

16( 



34 



2(M + 8) 2M 2 + 33M + 130 

17 ~ M+14 ' l3 ~ M2 + 22M + 112 (25j 

The two sets of expressions eqs. ( 1241(251) both come from identical Feynman diagrams. Note 
that the v exponents of 0$ and Sp^ agree with the identification iV = —M/2. This is 
to be expected since to this lowest order this exponent only involves the 1-point function of 
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XX an d doesn't involve a 2-point function of x-fields. The exponents (3 on the other hand 
do not agree under this identification. 

In the table below we summarize some of the above formulas and compare them with 



(3) 

the highest order epsilon expansions available to date. The column Q \J h igh — eps." is the 
Borel-summed 7-loop epsilon expansion for the 0$ model, taken from 
low — eps." is the lowest order epsilon expansion eq. (1251) . 
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24j. The column 



(3) 
M 




C„(3) 


Ojj high-eps. 


Of) low-eps. 


M = 


V 


8/7= 1.143 
33/28 = 1.178 


1.176 ± .002 
1.210 ± .003 


8/7= 1.143 
65/56 = 1.161 


3D Ising 
M = 1 


V 


5/4 = 1.250 
13/10 = 1.300 


1.261 ± .002 
1.303 ± .006 


6/5 = 1.200 
11/9 = 1.222 


M = 2 


V 


4/3 = 1.333 
25/18 = 1.388 


1.341 ± .003 
1.388 ± .006 


5/4 = 1.250 
51/40 = 1.275 


M = 3 


V 


7/5 = 1.400 
51/35 = 1.457 


1.415 ± .007 
1.465 ± .01 


22/7 = 1.294 
247/187= 1.321 


M = 4 


V 

(3 


16/11 = 1.454 
133/88 = 1.511 


1.482 ± .012 
1.532 ± .018 


4/3 = 1.333 
49/36 = 1.361 



The agreement of the first two columns is rather striking, but still, a few comments are 
in order. Most of our (3 exponents are within error bars. For the Ising case, (3 is within 
error bars and accurate to 1 part in 1000. The v exponents, though very reasonable, are 
not as accurate, and this may be due to the fact that we only calculated it to 1-loop, 
whereas j3 required a 2-loop calculation. (The RG beta-function is still only at 1-loop, 
however the (3 exponent depends on r y x and the lowest contribution to the latter arises at 2 
loops[3].) For the Ising model, much more is known, and with smaller errors, since it has 
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19 



20 



2 II ]. These works obtained 



been studied in a variety of Monte-Carlo simulations 
v = 1.2596, 1.2612, 1.2602, 1.2588 respectively, all with errors of about ±(.001 - .002). A 



non-perturbative RG approach gives v = 1.264 22j. So our 1-loop calculation of v is off by 



about 1% in the Ising case. We emphasize that this does not yet spoil our conjecture on the 
equivalence M = Sp_ 2M) rather it suggests that one must go to higher loop order to be 
within the error bars of what is already known. 
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IV. CONCLUSIONS 



The excellent agreement of the exponents of the N and Sp_2 N m °dels certainly strongly 
supports our conjecture. We also gave some arguments, though incomplete, towards an 
explanation based on symmetries, in particular the fact that O(N) is a subgroup of Sp(2N). 
Because the v exponent is so well measured for the Ising model, it appears our exponent is 
off by 1% in that case, but we suggested that this could be due to the fact that we only went 
to 1-loop for v. The f3 exponent on the other hand, which involves the 2-loop calculation 
of 7 X , is within the small error bars. Our calculation was a quite simple 1- and 2-loop 
calculation, so, if our conjecture is correct, the symplectic fermion is at least a very efficient 
way to compute exponents compared to the delicate Borel summations one has to do in the 
very high order epsilon expansions in the bosonic description. 

Higher order calculations will certainly change the exponents, so one may worry that the 
agreement so far is too good. However we wish to point out that for the Om models, it is 
known that when one goes to higher orders, the agreement becomes worse, which is why one 
needs to Borel sum. Higher order calculations are currently in progress and will be reported 
in a future publication. 
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